Combined CI+MBPT calculations of energy levels and transition amplitudes in Be, 

Mg, Ca, and Sr 
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Configuration interaction (CI) calculations in atoms with two valence electrons, carried out in the 
y(N-2) jjaj-tree-Fock potential of the core, are corrected for core-valence interactions using many- 
body perturbation theory (MBPT). Two variants of the mixed CI-I-MBPT theory are described and 
applied to obtain energy levels and transition amplitudes for Be, Mg, Ca, and Sr. 

PACS numbers: Sl.lO.+z, 31.25.-v, 32.30.-r, 32.70.-m 



I. INTRODUCTION 

Although Be, Mg, Ca, and Sr atoms have been studied 
theoretically for many years and numerous calculations 
are available in the literature, energy levels of those di- 
valent atoms have been treated primarily with semiem- 
pirical methods and only a limited number of low-lying 
levels have been evaluated using ab-initio methods, which 
often do not provide sufficient precision or require exten- 
sive computer resources. Semiempirical methods, to their 
advantage, do not require significant computer resources 
and can be applied easily to a large number of levels; 
however, such theories have limited predictive power and 
accuracy. Although energies obtained using semiempiri- 
cal methods agree well with one another and with exper- 
iment, oscillator strengths obtained by different semiem- 
pirical calculations are inconsistent Examples of 
semiempirical calculations can be found for Be in Q , for 
Ca in Ql, and for Sr in 0]. Large-scale ab-initio con- 
figuration interaction (CI) calculations of energies and 
transition rates, although capable of high accuracy, have 
been performed only for a few low-lying levels in the Be 
[Q, H and Mg |^ isoelectronic sequences. The size of the 
configuration space in such CI calculations is limited by 
the available computer resources. Smaller-scale CI cal- 
culations, carried out in the frozen Hartree-Fock 
potential of the core, lead to poor results. We found, for 
example, that frozen-core CI calculations in Ca gave en- 
ergies so inaccurate that it was difficult, if at all possible, 
to identify many closely spaced levels of experimental 
interest. Multi-configuration Dirac-Fock (MCDF) and 
Hartree-Fock (MCHF) methods have also been used to 
obtain energies and oscillator strengths in divalent atoms: 
MCHF for Be-like ions [7| and neutral calcium j^, and 
MCDF for Mg-like ions ||. The accuracy of MCHF and 
MCDF calculations in neutral atoms is poor, basically 
because of computational limits on the number of con- 
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figurations. Polarization potentials have been used in 
conjunction with MCHF calculations to improve the 
accuracy of energies for Ca I and Ca II. Many-body per- 
turbation theory (MBPT) calculations of energies and 
oscillator strengths for neutral divalent atoms using an ef- 
fective Hamiltonian within a small model space, are also 
found to be inaccurate |l^. Good agreement with 
experiment for divalent atoms, however, was achieved in 
Refs. m with a combined CI-hMBPT method. 

A related method was applied to calculations of ener- 
gies and oscillator strengths for Mg-like ions in Ref. [p^ . 
Among the ab-initio methods, CI-I-MBPT is particularly 
attractive since it is capable of giving accurate energies 
and transition rates for both light and heavy divalent 
atoms with modest computer resources. 

A precise and efficient theoretical method for cal- 
culations of properties of divalent atoms is needed 
for many possible applications of current interest, in- 
cluding calculations of spectra, transition amplitudes, 
hyperfine structure constants, polarizabilities, parity- 
nonconserving (PNC) amplitudes, van der Waals coef- 
ficients, and Lennard- Jones coefficients. There is also 
growing interest in properties of divalent atoms in con- 
junction with low-temperature Bose-Einstein condensa- 
tion (BEC) experiments. For example, the prospect for 
achieving BEC in divalent atoms was discussed in |jl^, |lj] 
and depends on the size of the van der Waals coefficient. 

At least two major difficulties have been recognized in 
studying divalent atoms. First, core polarization effects 
are significant and must be taken into account. A sim- 
ilar situation exists in monovalent atoms where various 
methods have been successfully applied to describe the 
valence-core interaction. We have made extensive use of 
one of these methods, MBPT, and have developed meth- 
ods for calculating all diagrams up to the third order 
for energies and transition amplitudes ||2^. A sec- 
ond major difficulty is that two valence electrons interact 
so strongly in neutral atoms that two-particle diagrams 
must be included to infinite order. Since infinite order is 
required, the MBPT method is difficult to apply. How- 
ever, valence-valence correlations can be accounted for 
completely using the CI method. 

With this in mind, we have developed a method (sim- 
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ilar to that used in Refs. |]T^ |T5j but with important 
differences) for high-precision calculations of properties 
of atoms with two valence electrons. The method starts 
with a complete CI calculation of the interactions be- 
tween the valence electrons in a frozen core and accounts 
for valence-core interactions using MBPT. We apply this 
combined CI-I-MBPT method to calculate energy levels 
and transition amplitudes for Be, Mg, Ca, and Sr. 



B. Combining CI with MBPT 

Core polarization effects can be treated using MBPT. 
In this paper, we introduce two procedures that enable us 
to combine frozen-core CI and second-order two- valence- 
electron MBPT, which we refer to as "CI averaging" and 
"Brueckner-Orbital CI" methods. 



CI averaging 



II. METHOD 



Frozen-Core CI 



We start with a lowest-order description of a divalent 
atom in which the closed N-2 electron core is described 
in the HF approximation and valence or excited elec- 
trons satisfy HF equations in the "frozen" T/P^"^) HF 
core. As we mentioned in the introduction, the strong 
valence-valence correlations must be included to infinite 
order; the CI method accomplishes this. The configura- 
tion space for divalent atoms is built up in terms of the 
excited HF orbitals. We include all orbitals with angular 
momentum / < 5 (partial wave contributions scale as 
+ 1/2)'*) and we use 25 basis functions out of a com- 
plete set of 40 for each value of angular momentum. The 
effect of these restrictions is insignificant considering the 
perturbative treatment of valence-core correlations. 

A detailed discussion of the CI method (as used here) 
can be found in Ref . |^ . We introduce a configuration- 
state wave function <!>/ = ^jM{ij) in which single- 
particle basis orbitals i and j are combined to give a 
two-particle wave function with angular momentum J 
and definite parity. We then expand the general two- 
particle wave function ^'jm in terms of all ^jAiiij) in 
our basis set 



(1) 



The expectation value of the Hamiltonian becomes 

(*JM Iff I *.7Af) = ^ C? + ^ ViK CI CK, (2) 



I,K 



where Ej = ei + ej is the sum of single-particle HF en- 
ergies and ViK is a first-order, t wo-p article correlation 
matrix element (see, for example, |2l)) between the con- 
figurations / = (ij) and K ~ [kl). The variational con- 
dition leads to CI equations 



^ {EjSiK + Vik) ck = Ac/, 



(3) 



K 



from which CI energies (A) and wave functions {J^i ci^i) 
are found. 



In this first method, the core- valence interaction AEyc 
is obtained by "averaging" MBPT corrections over CI 
wave functions: 



K 



(4) 



where the configuration weights c/ and ck are taken from 
the solution of the CI equation, Eq. (||), and i?'^-* is that 
part of the effective Hamiltonian projected onto the va- 
lence electron subspace containing second-order valence- 
core interactions. The dominant second-order parts of 
the effective Hamiltonian, beyond those accounted for in 
the CI calculation, are the screening and self-energy dia- 
grams: 77^^' = the self-energy being much 
larger than the screening and both being larger than the 
remaining second-order terms. 

We borrow ready-to-use formulas, derived using stan- 
dard techniques, from Ref. ||l^. The screening contribu- 
tion to the effective Hamiltonian is 



rrscrccn 
v' w' vw 



Ci{a'l3'al3) x 



E 

nbk 



[k] 



ja' jfS- J 
313 ja k 



Zk{a'ban)Zk{(3'nl3b) 



where 



Ci{a'P'aP) = (-1)"^ [Sa'v'S/S'w'SavSpw + Sa'w'S/3'^ 

The self-energy contribution to H'^^^ is 



TTSC 



self 

v' w'vw 



J 



(5) 



(6) 



(7) 



where 
Sij(eo) 



E 

kcmn 



(-1) 



X]^ (icmn) Z^. (jnnjc) 



^ (_ 1 ) J, -Je Xk {mbc) Zfc {bcjn) 



kbc 



Mk] 



eo + e„ - fib 



In the above equations, J is the angular momentum of 
the coupled two-particle states. The coupled radial inte- 
grals Xk{abcd) and Zk{abcd) are defined in ||l^. We use 
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the notation [k] = 2k + I. The quantities r]yw are nor- 
mahzation constants, r^i,^, = l/'v/Z for identical particle 
states and 1, otherwise. In the expression for the self- 
energy, the angular momenta of the ith and jth orbitals 
satisfy Ki = kj, where Ki = ^{ji -1-1/2) for ji = li zL 1/2 
is the angular quantum number uniquely specifying the 
spinor for state i. Since we found that the second-order 
self-energy correction is very important, we also consider 
the fourth-order self-energy obtained by iteration: 

S.,(eo)-S.,M+E^^^^^^^^^^^- (9) 

t^^ - 

In heavy atoms, the choice of eg deserves special consid- 
eration. Problems with denominators arise from the fact 
that single-particle orbitals used in the self-energy cal- 
culation are not optimal, in the sense that there is mu- 
tual interaction between valence electrons not accounted 
for, even approximately, in the potential and ac- 

counted for excessively in the V^^^ potential which is 
used, for example, in Ref. ||lj]. One practical solution to 
this problem is to use "optimized" denominators [|l4| . A 
consistent theory requires an ab-initio treatment of the 
denominator problem. Basing calculations of atoms with 
two valence electrons on a more realistic potential can 
reduce uncertainties in the choice of the denominator in 
the self-energy corrections. 

We calculated energies of several levels using the CI 
averaging method and found that the best agreement 
with experiment for Be and Mg was obtained with cq 
equal to 1/2 of the CI energy. For the case of Ca, the 
best agreement was obtained choosing eo between 1/2 
and 1 times the CI energy. One advantage of the CI av- 
eraging method is that the basic CI code is simple and 
that the CI wave functions can be stored and used many 
times. A cut-off condition can be imposed, as a compro- 
mise between speed and accuracy. The fastest approxi- 
mation (giving the poorest accuracy) is obtained by re- 
stricting the MBPT corrections to the leading configura- 
tions. We used this leading configuration approximation 
to estimate the magnitude of the core-excitation effects 
as the first step in developing our computer code. Ad- 
justing the cut-off condition, we readily reached a high 
level of accuracy (finally we chose the cut-off condition 
\ciCk\ < 0.002 for all calculations). The energies for 
several states of Be, Mg, and Ca presented in this paper 
have been calculated with the CI averaging method. The 
principal drawback of this method is that wave functions 
necessary for calculations of other properties are not au- 
tomatically obtained. 



2. Brueckner- Orbital CI 

The effective Hamiltonian formalism |Q leads to the 
problem of diagonalizing the Hamiltonian matrix built on 
the frozen-core two-electron configuration state functions 



$/. We split this matrix into functionally distinct pieces: 

where iJ*^"^ is the zeroth-order Dirac-Fock Hamiltonian, 
which in the Dirac-Hartree-Fock (DHF) basis is 

"-v'w'vw ~ "vv"Jiuw' y^v ~r J , 

and if*^^-' is the first-order electron-electron interaction 
Hamiltonian 

v'w'vw v'w'vw^ 

defined in Ref. p2| . iJ*^^^ is the second-order correction 
which consists of the two-particle screening correction 
and the one-particle self-energy correction defined pre- 
viously. In the CI averaging method, -hij(i) is diag- 
onalized first in a DHF basis (where TJ'"' is diagonal) to 
give state energies and CI wave functions, then H^"^^ is 
evaluated using the CI wave functions to give corrections 
for the core- valence interaction. 

In the Brueckner-orbital (BO) CI method, the ba- 
sis functions are chosen as orthonormal solutions of the 
quasi-particle equation, 

[/lo + V^HF + Sy(e)]0j = e^O0,. (11) 

In this BO basis, 

/^(O) ^ ^sclf\ ^ 5,,,^_,(,f o + ,^0). (12) 

\ / v'w'vw 

The basis orbitals include second-order self-energy cor- 
rections together with the lowest-order DHF potential. 
The residual non-trivial part of the effective Hamilto- 
nian in the BO basis is the sum H^^^ -\- iJ^^een ^j^g 
Bruckner orbital-CI method, the residual Hamiltonian 
matrix is evaluated in the BO basis and diagonalized to 
obtain state energies and CI wave functions. The BO- 
CI method is equivalent to CI averaging method if we 
neglect energy differences in the denominators of i^^^^f 
and i/scrccn Qj.(igj. Qf ij^g valence-valence interaction 
energy) , which are small compared to the core excitation 
energies. The BO-CI method is also equivalent to the 
effective Hamiltonian method in |Q to the same level 
of precision, provided all second-order diagrams are in- 
cluded. Some advantage is gained in accuracy compared 
to the CI averaging method, since the largest valence-core 
corrections [those from I]ij(eo)] are taken into account to 
infinite order. 

The Brueckner-orbital CI method is very convenient 
for calculations of transition amplitudes; once the resid- 
ual interaction is diagonalized, the associated wave func- 
tions are immediately available. We include random- 
phase approximation (RPA) corrections in calculations 
of transition amplitudes by replacing "bare" matrix el- 
ements with "dressed" elements as explained in [ po[ . 
Length-form and velocity-form dipole matrix elements 
are found to be in close agreement in BO-CI calculations 
that include RPA corrections. 



4 



TABLE I: Comparison of Cl-averaging energy levels (cm ^) 
of Be I with experimental data from the NIST database |p2[ . 



Config. 


Term 


J 


NIST 


Cl-average 


2s' 


IS 











2s2p 


3P° 





21978 


21996 


2s2p 


3P° 


2 


21981 


22000 


2s3s 


3S 


1 


52081 


52074 


2p2 


ID 


2 


56882 


56890 


2s3p 


3P° 


1 


58907 


58890 


2s3p 


3P° 


2 


58908 


58896 


2p2 


3P 


1 


59695 


59749 


2p2 


3P 


2 


59697 


59747 


2s3d 


3D 


3 


62054 


62033 


2s3d 


ID 


2 


64428 


64414 


2s4s 


3S 


1 


64506 


64528 


2s4s 


IS 





65245 


65261 


2s4p 


3P° 


2 


66812 


66792 


2s4d 


3D 


3 


67942 


67924 


2s4f 


3F° 


3 


68241 


68224 


2s4f 


1F° 


3 


68241 


68224 


2s4d 


ID 


2 


68781 


68774 


2s5s 


3S 


1 


69010 


69056 



III. CALCULATIONS OF SPECTRA USING CI 
AVERAGING 



TABLE 11: Comparison of frozen-core CI energies (cm ^) and 
Cl-averaging energies for Be I with experimental energies from 
the NIST database Q. 



Config 


Term 


J 


NIST 


Cl-average 


Diff. 


Frozen CI 


Diff. 


2s3s 


IS 





54677 


54664 


-13 


54509 


168 


2p2 


3P 





59694 


59737 


43 


60090 


-396 


2s5s 


IS 





69322 


69307 


-15 


69387 


+65 



A comparison of the resulting CI energies with mea- 
sured energies from the National Institute of Standards 
and Technology (NIST) database H is shown in Table |. 
This comparison provides the first test of the CI averag- 
ing method. The values listed in the table agree with 
experiment at the level of tens of cm^^. The residual 
deviation can be explained as neglect of small Coulomb 
and Breit diagrams, which will be the subject of future 
investigations. 

It is also interesting to compare CI energies, with and 
without the MBPT corrections AE^c, with energies from 
the NIST database. Such a comparison is given in Ta- 
ble H and illustrates the importance of the valence-core 
corrections. 

The agreement with experiment improves by an order 
of magnitude for the Cl-averaging method as compared 
with a frozen-core CI calculation. Indeed, we found it 
necessary to use the more precise energies obtained from 
the Cl-averaging method to properly identify the transi- 
tions shown in this table. 



The CI averaging method is fast and convenient for 
calculations of energies when a large number of levels are 
needed, especially at the stage of adjusting the code pa- 
rameters. Below, we present our calculations for many 
levels of Be, Mg, and Ca atoms to demonstrate the accu- 
racy of this method. We evaluate the valence-core cor- 
rection AEyc to the CI energy using a subset of the CI 
coefficients limited to those satisfying |c/Cif| < 0.002. 
The parameter eo in the self-energy was chosen to be 
eci/2 for Be and Mg. For calcium it was increased to 
3eci/4 to obtain better agreement for energies of the 4p^ 
states. 

The basis set used to set up the calculations consisted 
of 25/40 DHF basis functions for each value of ^ < 5. 
The basis functions were formed as linear combinations 
of B-splines of order 7, constrained to a cavity of radius 
i?=80 ao. 



A. Calculations for Be 

We chose to study a Be atom for several reasons. First, 
this atom has a small core and, consequently, requires 
relatively little computation time. Second, because of the 
small size of the core- valence interaction, calculations for 
Be are expected to be very precise. 



B. Calculations for Mg 

Another example where the CI averaging method pre- 
dicts energy levels accurately is magnesium. In this atom, 
however, core correlations are larger and the treatment 
of the valence-core interaction term requires more careful 
analysis. One important aspect is choosing the parame- 
ter eo in the denominators of the MBPT corrections, an- 
other is the treatment of self-energy diagrams. We found 
mild sensitivity of final energies in Mg to the choice of 
eo- The corrected energies shown in the column headed 
'CI -I- 2nd' in Table [II, which were obtained with the 
choice eo = eci/2, are seen to be in close agreement with 
experimental energies . 

Typically, the self-energy correction is much larger 
than other valence-core diagrams; for example, in the 
Mg ground state, the self-energy is —1.65 x 10~^ a.u. 
while the screening contribution is ten times smaller, 
1.83 X 10~^ a.u. Valence-core contributions in fourth- 
order, obtained by iterating (or chaining) the second- 
order Brueckner corrections are also found to be signifi- 
cant, —6.57 X 10~^ a.u. for the Mg ground state. The ef- 
fect of including corrections from chaining the self-energy 
shown in the column headed 'CI + 4th' in Table III is 



seen to further improve the agreement with experiment. 
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TABLE III: Comparison of energies (a.u.) in Mg obtained 
from frozen-core CI, Cl-averaging with 2nd-order self-energy, 
and Cl-averaging with chained 4th-order self-energy, with ex- 
perimental energies from the NIST database P2|. 



Conf. 


Level 


CI 


CI-h2nd 


CI-h4th 


Expt. 


A {cm-') 




'So 


0.818 


0.8329 


0.833513 


0.833518 


1 


3s4s 


'So 


0.624 


0.6349 


0.635260 


0.635303 


9 


3s5s 


\So 


0.583 


0.5938 


0.594240 


0.594056 


40 


3s6s 


'So 


0.566 


0.5772 


0.577813 


0.577513 


66 




'Po 


0.562 


0.5695 


0.569747 


0.570105 


79 


Ss3p 


'Pi 


0.723 


0.7336 


0.733991 


0.733869 


27 


3s3p 


'Pi 


0.661 


0.6733 


0.673673 


0.673813 


31 


3s4p 


'Pi 


0.604 


0.6156 


0.615834 


0.651524 


68 


3s4p 


'Pi 


0.597 


0.6086 


0.608606 


0.608679 


16 


3sSp 


'P2 


0.723 


0.7333 


0.733867 


0.733684 


67 








C. 


Ca atom 







In Table IV, several even parity J = levels are calcu- 
lated with the frozen-core CI and Cl-averaging methods. 
Compared to the frozen-core CI method, the agreement 
is significantly improved with the addition of MBPT cor- 
rections, changing the difference between experiment and 
theory from approximately one thousand cm^^ to a few 
hundred cm^^. This significant change clearly indicates 
the importance of the valence-core interaction, which is 
much stronger than in the case of Be and Mg. As a 
result, the final accuracy of CI-I-MBPT method is also 
lower than for the lighter atoms. While the poor accu- 
racy of frozen CI energies prevents the identification of 
energy levels, more accurate CI-I-MBPT energies permit 
one to identify many Ca levels. It is interesting to notice 
that the sequence of experimental levels for the states of 
a particular symmetry is the same as the sequence of the- 
oretical eigenvalues. Once the question of classification 
is solved, various properties of atoms can be calculated 
using, for example, frozen-core CI. 

In the case of Ca, another problem that needs atten- 
tion is the choice of the parameter eo in the self-energy. 



TABLE IV: Comparison of the accuracy of frozen-core CI 
and CI averaging calculations for Ca. The parameter eo — 
0.75 eg. 

Conf. Level frozen CI Diff. Cl-average Diff. Expt. 



4s5s 


'So 


31901 


-1416 


33196 


-121 


33317 




'Po 


36699 


-1718 


38900 


483 


38418 


4s6s 


'So 


39376 


-1314 


40504 


-186 


40690 


4p2 


'So 


41480 


-306 


42366 


580 


41786 


4s7s 


'So 


42673 


-1604 


43841 


-436 


44277 


4s8s 


'So 


44277 


-1610 


45551 


-336 


45887 


4s9s 


'So 


45629 


-1206 


46912 


77 


46835 



TABLE V: Comparison of DHF spline energies 'DHF', 
second-order energies '2nd order', and energies resulting from 
diagonalization of the self-energy matrix, Brueckner-orbital 
'BO' energies, with experiment for the Mg+ ion. The core 
configuration is ls^2s^2p^ . The size of the self-energy ma- 
trix is 14x14 for each angular momentum. All energies are 
expressed in cm~^. 



States 


DHF 


2nd order 


BO 


Expt. 


3Si/2 


118825 


121127 


121184 


121268 


4Si/2 


50858 


51439 


51446 


51463 


5Si/2 


28233 


28467 


28469 


28477 


3pi/2 


84295 


85508 


85542 


85598 


4pi/2 


40250 


40625 


40633 


40648 


5pi/2 


23642 


23808 


23811 


23812 



the dominant part of the core-valence interaction. We 
find that there is an optimal value of this parameter be- 
tween ec7/2, our standard value for Be and Mg, and 
eci, for which the ground state becomes very accurate. 
In Table |^ we chose this parameter to be 0.75 ec/. In 
the following section, we will illustrate our calculations of 
transition amplitudes for several levels of Mg, Ca, and Sr 
where other precise calculations and measurements exist. 



IV. CALCULATIONS USING THE 
BRUECKNER-ORBITAL CI METHOD 

In this section, we present our calculations of energies 
and transition amplitudes with the Brueckner-orbital CI 
method. Our basis consisted of 25 V^f^ orbitals (those 
orbitals were constructed of 40 B-splines in the cavity 
80 a.u.), in which 14 lowest excited states were replaced 
with Brueckner orbitals. The resulting one-valence elec- 
tron energies for the divalent atoms were tested by com- 
paring with experimental energies for the corresponding 
monovalent ions. For Mg+, the BO energies agree with 
experiment better than do the second-order energies (Ta- 
ble ^ . A second iteration of the BO equation was also 
included in the Cl-averaging method (Table III) to im- 
prove accuracy. The small size of the residual devia- 
tion from experiment in both tables can be attributed 
to higher-order diagrams. Two-particle screening cor- 
rections with the restriction n < 15 were included in 
the effective Hamiltonian, diagonalization of which pro- 
vided the initial and final state wave functions necessary 
for the calculation of transition amplitudes. We checked 
that restrictions on the number of BO and screening di- 
agrams included in the calculation did not lead to sig- 
nificant errors. Dressed transition amplitudes were used 
to take into account RPA corrections, which provide bet- 
ter length- and velocity- form agreement. We completely 
neglected the extremely time consuming structural ra- 
diation corrections which are expected to be small for 
the length form; for this reason, the result calculated 
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TABLE VI: Comparison of the present transition energies uj 
(a.u.) and oscillator strengths / for Be with those from other 
theories and experiment. A few allowed singlet-singlet transi- 
tions of the type So — P° between low-lying states are consid- 
ered. The experimental uncertainties are given in parentheses. 



TABLE VII: Comparison of BO-CI energies (cm ^) with ex- 
periment for Mg, Ca, and Sr. 



Levels 



Theory 



Expt. 



Diff. 



Transition 


Source 


(Theory) 


a;(Expt.) 


f 


2s^-2s2d 


present 


0.194126 


193954 


1.3750 






24 








1 38(0 121 






25 








1 34C0 05"! 






iL 


0.91412 




1.375 






2J 




193914 




1.374 






26 








1.3847 






27 








1.470 






~ 








1.375 


2s^-2s3p 


present 


0.274231 


0.274251 


00904 








0.27441 




n 00Q01 






2o 


1 


0.274236 




0.00914 














0.0104 






27 








0.037 












n 00885 


2s3s-2s2p 


present 


0.054977 


0.05519 


0.1188 








0.05509 




0.118 






23 




0.055198 




0.1175 






26 


1 






0.1199 






27 








0.140 


2s3s-2s3p 


present 


0.025128 


0.025107 


0.9557 






24 




0.0252 




0.958 










0.025124 




0.9565 






26 








0.9615 



in length form should be considered as more accurate. 
Small normalization corrections are also omitted. 



A. Be case 

The most accurate results for divalent atoms are ex- 
pected for Be since it contains the smallest MBPT cor- 
rections. In Table Vl, we compare our calculations with 
available precise calculations and experiment. Transi- 
tion energies agree with experiment to better than 0.1%, 
except for the transition 23^8^3 — 2s2p^P which has 
0.4% accuracy. Our oscillator strengths agree well with 
those obtained in very accurate ab-initio calculations of 
Ref. 1^ and in semiempirical calculations of Ref. jl| that 
reproduce energies very closely; for the principal transi- 
tion 2s2^S - 2s2p^P, our value 1.375 differs by 1 in the 
4th digit from the value 1.374 in Ref. the accuracy 
being better than 0.1%, and coincides with the value of 
Ref. [|l|. Very close agreement with ab-initio theory is 
also achieved for the transition 2s3s^S — 2s3p^P. For 
suppressed transitions, an accuracy of 1% is obtained. 
Conducting a simple statistical analysis, we found that 



3s4s ^So 
3s5s ^So 
3s6s ^So 
3s3p ^'Pi 
3s3p ^Pi 
3s4p ^'Pi 
3s4p ^Pi 

4s5s ^So 
4p2 'Po 
4s6s ^5*0 
4s4p ^Pi 
4s4p ^Pi 
4s5p ^Pi 
4s5p ^Pi 

5s6s ^So 
5p' 'Po 
5p' iPo 
5s5p ^Pi 
5s5p ^Pi 
5s6p ^Pi 
5s6p ^Pi 



Mg atom 
43452 
52517 
56154 
21834 
35059 
47806 
49317 

Ca atom 
33505 
38651 
40862 
15595 
23797 
36760 
36917 

Sr atom 
30874 
35913 
37696 
15081 
21981 
34293 
34512 



43503 
52556 
56187 
21870 
35051 
47844 
49347 

33317 
38418 
40690 
15210 
23652 
36555 
36732 

30592 
35193 
37160 
14504 
21699 
33868 
34098 



-51 
-39 
-33 
-44 
8 
-38 
-30 

188 
233 
172 
385 
145 
205 
185 

282 
720 
536 
577 
282 
425 
414 



energy differences in the Cl-averaging and BO-CI calcula- 
tions have similar statistical errors, but slightly different 
systematic shifts which can be explained partially by dif- 
ferent denominators in the two methods. Another reason 
is the cut-off condition 0.002 in the former method and 
restriction on the number of Brueckner orbitals in the 
latter. The effect of the partial wave restriction on the 
ground state energy in both methods is 6 cm~^. If this 
value is accounted for, the agreement becomes slightly 
better. The results in our tables are not extrapolated 
owing to the smallness of the omitted partial wave con- 
tributions. 



B. The cases of Mg, Ca, and Sr 

The accuracy of both the Cl-averaging and the BO- 
CI calculations considered above decreases from light to 
heavy divalent atoms. Table VII illustrates this tendency 
in BO-CI calculations: for Mg, the theory-experiment 
differences ra nge within 50 cm^^, similar to what we 



have in Table III, and for Ca the deviation from exper- 
iment increases to about 200 cm~^ which is comparable 



to that in Table IV. The lowest accuracy is for Sr, which 
has the largest core and MBPT corrections. Similar re- 
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suits for energies have been obtained in Ref.||T^. Our 
experiment-theory diflerences exhibit a systematic shift, 
which if subtracted, brings results into better agreement. 
For example, in Ca this shift is 216 cm"'^. After its sub- 
traction, the residual deviation is 73 cm^^. This subtrac- 
tion procedure can be used in cases where closely spaced 
levels are difficult to identify. The systematic shift can 
be attributed to omitted correlations that affect mostly 
the ground state which is used as a reference. The cut-off 
condition in the Cl-averaging method and restrictions on 
the number of BO and screening diagrams also has some 
effect on the accuracy of our results. This is one reason 
why the two methods give slightly different energies. In 
future development of our computer code, we will try to 
remove such restrictions completely. Another reason why 
the two methods give differe nt re sults is that the choices 
of £0 were different. In Table VIII , we illustrate our calcu- 
lations of transition amplitudes for Mg, Ca, Sr. All of our 
transition amplitudes completely agree with those of re- 
cent CI+MBPT calculations by Porsev et al. , and are 
close to experimental values. Length-form and velocity- 
form amplitudes agree to better than 1% for allowed tran- 
sitions. Forbidden transitions are more problematic, ow- 
ing to cancellation effects, and have poorer agreement 
between gauges and with experiment. The inclusion of 
the Breit interaction and negative-energy contributions, 
which are more important for the velocity form, might 
improve the situation. We also noticed that, if the bal- 
ance between states such as Pi/2 and P3/2 in relativistic 
basis is not properly maintained, the results for nonrela- 
tivistically forbidden transitions will be unstable. In ad- 
dition, those transitions were affected by the number of 
BO and screening diagrams included in calculations. To 
minimize or exclude those effects in the BO-CI method, 
the BO orbitals and cut-off conditions were made com- 
pletely symmetric with respect to Z + 1/2 and 1 — 1/2 
orbitals and included BO and screening corrections with 
number of excited orbitals less than 15. 



exist. Further improvement in accuracy is possible and 
is being pursued. This theory can be extended easily 
to treat particle-hole excited states of closed-shell atoms, 
atoms with three valence electrons, and other more com- 

TABLE VIII: Comparison of our length-form (L) and 
velocity-form (V) calculations with those from Ref. hm and 
with experiment. 





Mg 


Ca 


Sr 




1 no / \ 

P"{nsnp) 






L 


4.026 


4.892 


5.238 


V 


4.019 


4.851 


5.212 


Other" 


4.03(2) 


4.91(7) 


5.28(9) 


Expt. 


4.15(10)* 


4.967(9)' 


5.57(6)^ 




4.06(10)" 


4.99(4)^ 


5.40(8)'' 




4.12(6)'' 


4.93(11)» 












L 


0.0063 


0.0323 


0.164 


V 


0.0070 


0.0334 


0.166 


Other" 


0.0064(7) 


0.034(4) 


0.160(15) 


Expt. 


0.0053(3)' 


0.0357(4)' 


0.1555(16)° 




0.0056(4)^ 


0.0352(10)™ 


0.1510(18)™ 




0.0061(10)'= 


0.0357(16)" 


0.1486(17)" 



"Porsev et al. ||l 4] . 
''Liljcby et al. [^3. 
■^Lundin et al. | pC| |. 
''Smith and Gallagher H. 
■^Zinner et al. jis) . 
•'Kelly and_Mathur |32|. 
9Hansen [gj. 
''Parkinson et al. Q . 
'Godone and Navero jssj. 
'Kwong et al. Q . 
'^Mitchell ||^. 
'Husain and Robert 
'"Drozdowski et al. | 
"Whitkop and Wiesenfi 
"Husain and Schifino [ |l[ 
PKelly et al. M. 



V. SUMMARY AND CONCLUSION 

In this paper, we have introduced two methods to im- 
prove the accuracy of the frozen-core CI calculations us- 
ing MBPT: the Cl-averaging method and the Brueckner- 
orbital CI method. We have applied these methods to Be, 
Mg, Ca, and Sr atoms. Our calculated energies and tran- 
sition amplitudes for those atoms are in close agreement 
with the results of the best available theories and experi- 
ments. Compared to semiempirical theories, our method 
has an advantage in accuracy, and compared to other 
ab-initio theories, an advantage of simplicity. These two 
methods can also be used to evaluate properties of Ry- 
dberg states for which only semiempirical calculations 



plicated systems. 



Acknowledgments 

The authors are grateful to U. I. Safronova for help- 
ing to establish the correctness of computer codes. We 
are thankful to M. Kozlov for discussion of the theory. 
We thank H. G. Berry for reading manuscript and giving 
useful comments. We thank A. Derevianko for pointing 
out about experimental interest in divalent atoms. This 
work was supported in part by National Science Founda- 
tion Grant No. PHY-99-70666. 



[1] M.-K. Chen, J. Phys. B 31, 4523 (1998). 



[2] J. Mitroy, J. Phys. B 26, 3703 (1993). 



8 



C. J. Dai, Phys. Rev. A 52, 4416 (1995). 
M. H. Chen and K. T. Cheng, Phys. Rev. A 55, 166 [24 
(1997). 

M. H. Chen, K. T. Cheng, and W. R. Johnson, Phys. [25 
Rev. A 64, 042507 (2001). 
M. H. Chen and K. T. Cheng, Phys. Rev. A 55, 3440 [26 

(1997) . [27 
C. F. Fischer, M. Godefroid, and J. Olsen, J. Phys. B 30, 

1163 (1997). [28; 
N. Vaeck, M. Godefroid, and J. E. Hansen, J. Phys. B 
26, 361 (1991). 

P. Jonsson and C. F. Fischer, J. Phys. B 30, 5861 (1997). [29 
T. Brage and C. F. Fischer, Physica Scripta 48, 533 

(1993). [30 
W. R. Johnson, M. S. Safronova, and U. I. Safronova, 
Physica Scripta 56, 252 (1997). [31 
M. S. Safronova, W. R. Johnson, and U. Safronova, Phys. 
Rev. A 53, 4036 (1996). [32 
V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, Phys. 
Rev. A 54, 3948 (1996). [33 
S. G. Porsev, M. G. Kozlov, Y. G. Rakhlina, and A. Dcre- [34 
vianko, Phys. Rev. A 64, 012508 (2001). 

V. A. Dzuba and W. R. Johnson, Phys. Rev. A 57, 2459 [35 

(1998) . [36^ 
R. K. Chaudhuri, B. P. Das, and K. F. Freed, J. Chem. 
Phys. 108, 2556 (1998). [37 
M. Machholm, P. S. JuUenne, and K.-A. Suominen, Phys. [38 
Rev. A 59, R4113 (1999). 
G. Zinncr, T. Binnewies, F. Riehle, and E. Tiemann, [39 
Phys. Rev. Lett. 85, 2292 (2000). 

S. A. Blundell, W. R. Johnson, and J. Sapirstein, Phys. [40 
Rev. A 42, 3751 (1990). 

I. M. Savukov and W. R. Johnson, Phys. Rev. A 62, [41 
052512 (2000). 

M. H. Chen, K. T. Cheng, and W. R. Johnson, Phys. [42 
Rev. A 47, 3692 (1993). 

Available at http://physics.nist.gov/cgi-bin/AtData. 
K. T. Chung and X.-W. Zhu, Phys. Rev. A 48, 1944 



(1993). 

S. Hontzeas, I. Martinson, P. Erinan, and R. Buchta, 
Physica Scripta 6, 55 (1972). 

I. Martinson, A. Gaupp, and L. J. Curtis, J. Phys. B 7, 
L463 (1974). 

R. Moccia and P. Spizzo, J. Phys. B 18, 3537 (1985). 

E. Markiewicz, R. P. McEachran, and M. Cohen, Physica 
Scripta 23, 828 (1981). 

J. Fleming, M. R. Godefroid, K. L. BeU, A. Hibbert, 
N. Vaeck, J. Olsen, P. Jonsson, and C. F. Fischer, J. 
Phys. B 29, 4347 (1996). 

L. Liljeby, A. Lindgard, and S. M. et al, Physica Scripta 
21, 805 (1980). 

L. Lundin, B. Engman, J. Hilke, and I. Martinson, Phys- 
ica Scripta 8, 274 (1973). 

W. W. Smith and A. Gallagher, Phys. Rev. 145, 26 
(1966). 

F. M. Kelly and M. S. Mathur, Can. J. Phys. 58, 1416 
(1980). 

W. Hansen, J. Phys. B 16, 2309 (1983). 

W. H. Parkinson, E. M. Reeves, and F. S. Tomkins, J. 

Phys. B 9, 157 (1976). 

A. Godone and C. Novero, Phys. Rev. A 45, 1717 (1992). 
H. S. Kwong, P. L. Smith, and W. H. Parkinson, Phys. 
Rev. A 25, 2629 (1982). 

C. J. Mitchell, J. Phys. B 8, 25 (1975). 

D. Husain and G. J. Roberts, J. Chem. Soc. Faraday 
Trans. 2 82, 1921 (1986). 

R. Drozdowski, M. Ignasiuk, J. Kwela, and J. Heldt, Z. 
Phys. D 41, 125 (1997). 

P. G. Whitkop and J. R. Wiesenfeld, Chem. Phys. Lett. 
69, 457 (1980). 

D. Husain and J. Schifino, J. Chem. Soc. Faraday Trans. 
2 80, 321 (1984). 

J. F. Kelly, M. Harris, and A. Gallagher, Phys. Rev. A 
37, 2354 (1988). 



